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We  have  analyzed  the  colorimetric  and  spectral  characteristics  of  2600  daylight  spectra  (global  spectral  irra- 
diances  on  a  horizontal  surface)  measured  for  all  sky  states  during  a  2-year  period  at  Granada,  Spain.  We 
describe  in  detail  the  chromaticity  coordinates,  correlated  color  temperatures  (CCT),  luminous  efficacies,  and 
relative  UV  and  IR  contents  of  Granada  daylight.  The  chromaticity  coordinates  of  Granada  daylight  lie  far 
above  the  CIE  locus  at  high  CCTs  (>9000  K),  and  a  CCT  of  5700  K  best  typifies  this  daylight.  Our  principal- 
components  analysis  shows  that  Granada  daylight  spectra  can  be  adequately  represented  by  using  six¬ 
dimensional  linear  models  in  the  visible,  whereas  seven-dimensional  models  are  required  if  we  include  the  UV 
or  near-IR.  Yet  on  average  only  three-dimensional  models  are  needed  to  reconstruct  spectra  that  are  colori- 
metrically  indistinguishable  from  the  original  spectra.  ©  2001  Optical  Society  of  America 
OCIS  codes:  010.1290,  330.1710,  330.1730. 


1.  INTRODUCTION 

Detailed  knowledge  of  the  spectral  power  distribution 
(SPD)  of  daylight  in  the  visible,  ultraviolet,  and  near- 
infrared  parts  of  the  spectrum  at  different  sites  and  for 
varying  atmospheric  conditions  has  many  applications. 
Such  knowledge  is  important  in  scientific  disciplines  (e.g., 
color  vision,  atmospheric  optics,  meteorology,  medicine, 
biology,  and  materials  research),  technical  areas  (e.g., 
photovoltaic  applications,  biomass  production,  color  ren¬ 
dering  and  metamerism,  agriculture,  and  architecture), 
and  in  industry  (e.g.,  photography,  dyes,  paints,  and  tex¬ 
tiles). 

During  the  1960’s  and  1970’s,  measurement  campaigns 
were  conducted  in  several  countries  (United  States,  En¬ 
gland,  Japan,  India,  South  Africa,  Australia)  in  efforts  to 
determine  representative  daylight  spectra  E(\).1-14  Re¬ 
searchers  variously  made  these  measurements  using 
spectroradiometers, 1-2,4-5,7-14  colorimeters,3  or  visual 
color-temperature  meters.6  Furthermore,  their  measure¬ 
ments  involved  different  detector  orientations,  fields  of 
view,  and  sky  states,  and  they  sometimes  analyzed  day¬ 
light  and  skylight  data  simultaneously.  Despite  this  va¬ 
riety  of  techniques  and  equipment,  these  studies  consis¬ 
tently  reached  two  basic  conclusions.  First,  daylight’s 
different  phases  have  chromaticities  that  lie  near  the 
Planckian  locus  of  the  1931  chromaticity  diagram 
adopted  by  the  Commission  Internationale  de  l’Eclairage 
(CIE).  Second,  these  different  daylight  power  spectra  are 
highly  correlated  with  one  another,  and  this  underlying 
similarity  has  many  practical  uses.  However,  none  of 
these  older  studies  include  daylight  spectra  measured 
over  both  a  wide  range  of  atmospheric  conditions  and  a 
long  time  period.  For  example,  Condit  and  Grum  mea- 
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sured  251  spectra  in  just  four  weeks,4  and  Tarrant  mea¬ 
sured  only  391  spectra  over  four  years.11 

The  paper  by  Judd  et  al.  ,2  which  analyzes  three  sets  of 
daylight  spectra,  is  the  one  most  often  cited  in  the  litera¬ 
ture  and  is  the  basis  for  the  CIE’s  daylight  recommenda¬ 
tions.  Notable  among  the  results  of  Judd  et  al.  is  that 
most  daylight  SPDs  can  be  accurately  estimated  by  lin¬ 
early  combining  three  fixed  basis  functions.  Although 
Judd  et  al.  remains  perhaps  the  best  experimental  analy¬ 
sis  of  daylight  to  date,  its  results  depend  heavily  on  the 
particular  conditions  under  which  its  spectra  were  ac¬ 
quired,  such  as  sky  state,  field  of  view,  and  detector  ori¬ 
entation. 

Other  papers  on  daylight  and  skylight  spectra  have  re¬ 
peatedly  shown,1-14  in  agreement  with  Judd  et  al.,  that 
chromaticity  provides  good  estimates  not  only  of  day¬ 
light’s  visible  spectra  but  also  of  its  correlated  color  tem¬ 
perature  (CCT).  Such  findings  prompted  the  CIE  to  pro¬ 
pose  a  method  for  calculating  the  relative  SPD  of  typical 
daylight  phases  that  uses  only  CCT  values.15  When  the 
spectral  region  is  broad  (including  the  near-IR,  where 
there  is  lower  spectrum-to-spectrum  correlation),  Slater 
and  Healey  show  that  highly  accurate  daylight  spectra  re¬ 
quire  linear  models  of  seven  dimensions  rather  than  the 
CIE’s  three.16  In  addition,  Romero  et  al.  use  seven¬ 
dimensional  models  to  accurately  reconstruct  SPDs  of 
both  daylight  and  artificial  illuminants.17 

We  have  earlier  noted  the  potentially  confusing  nature 
of  daylight  terminology,18  and  others  have  attempted  to 
eliminate  this  confusion.19  To  avoid  ambiguity,  one  must 
specify  what  is  being  measured:  the  natural-light  com¬ 
ponent  (skylight,  direct  sunlight,  or  their  combination), 
radiometric  variable  (radiance  or  irradiance),  and  any 
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other  parameters  which  affect  SPDs  (e.g.,  weather,  detec¬ 
tor  orientation,  field  of  view,  sampling  interval).  In  ad¬ 
dition,  we  follow  Henderson’s  advice  that  it  is  better  to 
analyze  “homogeneous  sets  of  observations  for  single  lo¬ 
cations  than  to  lose  the  small  differences  by  combining  a 
large,  heterogeneous  collection  of  curves”  (Ref.  20,  p.  184). 

Although  daylight  is  essentially  achromatic,  many  fac¬ 
tors  (e.g.,  solar  elevation,  site  altitude,  atmospheric  con¬ 
ditions,  pollution,  detector  orientation)  can  affect  its  spec¬ 
tral  composition  and  thus  its  color.  For  our  purposes, 
analyzing  only  the  daylight  spectra  themselves  runs  the 
risk  of  hiding  perceptually  interesting  facts,  and  this 
leads  us  to  colorimetric  analysis  of  our  spectra.  As  Hend¬ 
erson  recommends,  we  acquired  daylight  spectra  over  a 
long  time  period  at  one  site,  and  this  makes  our  colorim¬ 
etry  both  novel  and  especially  useful.  Our  analysis  below 
(Section  3)  considers  the  chromaticity,  CCT,  photometry, 
and  relative  UV  and  IR  contents  of  our  spectra. 

To  optimize  artificial-vision  algorithms  for  recognizing 
and  identifying  colors,  researchers  require  rigorous  math¬ 
ematical  analysis  of  the  daylight  SPDs  that  illuminate 
objects.21,22  Yet  before  our  work,  suitably  large  data¬ 
bases  of  daylight  spectra  were  not  readily  available.  Our 
collection  in  fact  lets  us  do  a  principal-value  decomposi¬ 
tion  (Section  4)  that  determines  the  dimensionality  of  ac¬ 
curate  linear  models  of  these  spectra. 

2.  EXPERIMENTAL  PROCEDURE 

Here  our  SPDs  are  of  hemispheric  daylight:  global  spec¬ 
tral  irradiances  E( A.)  on  a  horizontal  surface  from  direct 
sunlight  (when  present)  and  the  entire  sky.  Over  a  pe¬ 
riod  of  2  years,  we  recorded  2600  spectra  from  the  roof  of 
the  University  of  Granada’s  Science  Faculty  (37°11'N, 
3°37'W,  altitude  680  m),  which  is  in  an  urban,  nonindus¬ 
trial  area.  The  spectra  used  here  are  part  of  a  larger 
project23  and  have  been  used  in  a  variety  of  ways,  includ¬ 
ing  developing  an  improved  algorithm  for  calculating  day¬ 
light  and  skylight  CCTs. 24,25 

Our  observing  site’s  horizon  has  few  obstructions,  and 
so  our  spectra  sample  nearly  2  it  sr  of  the  sky.  Specifi¬ 
cally,  nearby  buildings  project  no  more  than  9°  above  the 
astronomical  horizon  and  obscure  only  ~7%  of  the  sky. 
Given  our  detector’s  cosine  response,  the  effect  on  our 
measurements  is  insignificant.  We  used  a  LI-COR  model 
LI- 1800  spectroradiometer26  equipped  with  a  holographic 
grating  monochromator  and  a  cosine-corrected  receptor  to 
measure  E(\)  from  300  to  1100  nm  in  5-nm  steps  (this 
spectroradiometer  is  described  in  more  detail  else¬ 
where27"29).  We  calibrated  the  LI-1800  every  six  months 
with  a  200-W  quartz  tungsten  halogen  lamp,  and  we 
avoided  any  strong  temperature  variations  that  could 
change  its  photodetector’s  sensitivity. 

From  February  1996  to  February  1998,  we  made  mea¬ 
surements  of  E(\)  in  nearly  all  weather  (clear  sky,  partly 
cloudy,  mostly  cloudy,  hazy,  overcast,  misty,  etc.),  exclud¬ 
ing  only  rainfall  and  snowfall.  We  acquired  data  from 
sunrise  to  sunset  at  1-h  intervals,  switching  to  10-min  in¬ 
tervals  when  solar-elevation  angles  h0  <  +5°.  Our 
measurements’  minimum  h0  ~  —4°,  because  at  this  and 
smaller  h0  the  LI- 1800  has  unacceptably  low  signal-to- 
noise  ratios.  Simultaneously,  we  used  a  pyranometer  to 


measure  how  stable  the  broadband  global  irradiance  was 
during  a  spectral  scan:  when  integrated  irradiance  var¬ 
ied  significantly  (>1%),  we  rejected  the  corresponding 
E (A.)  scan. 

3.  EXPERIMENTAL  MEASUREMENTS 

Figure  1  shows  three  different  daylight  spectra  measured 
in  Granada,  each  normalized  at  560  nm.  The  daylight 
phase  whose  CCT  =  5770  K  corresponds  to  a  clear  day 
with  h0  =  69°.  The  curve  with  CCT  =  4250  K  was  mea¬ 
sured  during  an  overcast  sunrise,  and  the  24,380-K  spec¬ 
trum  was  taken  during  a  partly  cloudy  sunrise. 
Throughout  our  analysis  we  calculate  CCTs  by  using  the 
binary  search  algorithm  described  in  Ref.  25. 

A.  Chromaticity  Coordinates 

Figure  2  shows  the  CIE  1931  chromaticities  of  all  2600 
Granada  daylight  spectra  overlaid  on  the  Planckian  locus 
and  the  CIE  daylight  locus.20  Also  plotted  in  Fig.  2  is  our 
Granada  daylight  locus,  which  we  calculate  from  a  least- 
squares  fit  to  be 

y  =  -0.05188  +  1.55320*  -  1.09234*2.  (1) 

Our  Granada  daylight  spectra  provide  a  measured 
chromaticity  gamut18  unprecedented  in  the  daylight 
literature.20  Most  Granada  data  with  low  CCTs  lie  above 
the  Planckian  locus  and  below  the  CIE  locus,  in  agree¬ 
ment  with  Judd  et  al.2  However,  the  two  daylight  loci 
differ  markedly:  Like  other  researchers,  we  find  that  our 
locus’s  slope  is  smaller  than  that  of  the  CIE 
locus.1,3"5,9,10,14  For  high  CCTs  (>9000  K)  the  CIE  locus 
does  not  adequately  represent  Granada  daylight  chroma¬ 
ticities,  which  are  perceptibly  greener  than  the  CIE’s. 
All  Granada  chromaticity  coordinates  on  the  green  side  of 
the  CIE  locus  (at  CCT  >  9000  K)  are  generated  by  clear 
skies  and  small  h  0 . 

The  Granada  locus,  like  the  locus  derived  by  Condit 
and  Grum,4  is  farthest  from  the  Planckian  locus  at  high 
CCTs  (i.e.,  toward  the  greens).  In  fact,  the  greatest  1931 
x,  y  Euclidean  distance  between  the  Granada  and  CIE 
daylight  loci  is  —0.011,  or  some  seven  just-noticeable  dif¬ 
ferences  (JNDs).30  Our  new  daylight  locus  describes  its 
associated  chromaticities  quite  well.  In  other  words,  our 
data  scatter  very  little  about  it  despite  the  length  of  our 


Fig.  1.  Relative  spectral  irradiance,  normalized  at  560  nm,  for 
three  different  daylight  SPDs  measured  in  Granada. 
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Fig.  2.  CIE  1931  chromaticities  of  our  2600  Granada  daylight  measurements  (circles)  overlaid  on  the  CIE  daylight  locus  (dashed  curve), 
the  Planckian  locus  (solid  curve  with  squares),  and  the  Granada  daylight  locus  (dotted  curve).  The  inset  shows  the  entire  CIE  1931 
diagram  and  Planckian  locus. 


measurement  campaign  and  the  number  of  irradiance 
spectra  measured.  This  small  scatter  about  a  daylight 
locus  is  comparable  to  that  found  by  some 
authors1,3"5’7'10,14  but  is  much  less  than  that  observed  by 
others.911"1'3  For  example,  Sastri  and  Das9  did  not  cal¬ 
culate  a  daylight  locus  because  their  chromaticities 
spanned  so  wide  an  area  on  the  CIE  diagram. 

Most  of  our  daylight  chromaticities  are  above  the 
Planckian  locus  toward  the  greens.  In  fact,  only  30  of  our 
2600  chromaticities  lie  below  it  toward  the  purples.  Al¬ 
though  this  distribution  has  often  been  seen 
before,1’3’7’1011  some  southern-hemisphere  researchers 
find  that  most  of  their  daylight  chromaticities  are  below 
the  Planckian  locus.9’13’14  To  our  knowledge,  no  one  has 
satisfactorily  explained  this  difference  between  hemi¬ 
spheres,  although  the  wide  variety  of  sites  in  both  sug¬ 
gests  that  local  terrain  color  is  not  responsible. 

B.  Correlated  Color  Temperatures 

As  noted  above,  researchers  have  consistently  shown  that 
daylight  CCTs  allow  accurate  estimates  of  its  visible- 
wavelength  irradiance  spectra  E(X).  This  makes  CCT  a 
useful  shorthand  for  specifying  daylight’s  colorimetric 
and  spectral  characteristics.  In  fact,  the  CIE  describes 
the  relative  spectral  irradiances  of  typical  daylight  phases 
at  the  Earth’s  surface  by  invoking  CCTs  ranging  from 
4000  to  25,000  K.15 

The  CCT  limits  for  our  2600  daylight  measurements 
are  3758-34,573  K.  To  make  comparing  our  data  with 
that  of  other  researchers  more  perceptually  meaningful, 
we  use  inverse  CCT  rather  than  CCT  itself.  The  unit  of 
inverse  CCT  is  the  reciprocal  mega-Kelvin  ( 106/CCT),  de¬ 
noted  by  the  symbol  MK"1  (historically  known  as  the 
“mired”).31  Reciprocal  mega-Kelvins  produce  a  uniform 


Fig.  3.  Histogram  of  inverse  correlated  color  temperature  (CCT) 
for  our  2600  Granada  daylight  spectra.  Each  bin  is  5  MK"1 
wide. 

scale  that  better  describes  human  color  sensitivity  than 
does  CCT  proper.  On  this  inverse  scale,  the  inverse 
CCTs  most  frequently  produced  by  our  Granada  daylight 
spectra  range  from  175  to  180  MK"1  (or  5555  to  5715  K), 
as  shown  in  Fig.  3. 

Table  1  compares  the  most  frequently  observed  (or 
modal)  intervals  of  inverse  CCT  from  several  researchers’ 
daylight  measurement  campaigns.1,6’711’14’20  Call  this 
modal  interval  m(CCT).  Only  clear-sky  hemispheric  ir¬ 
radiances  have  a  consistent  m(CCT),  and  it  is  —170-180 
MK"1.  Note  that  we  include  only  a  few  papers  in  Table  1 
because  of  (1)  intrinsic  differences  caused  by  differences 
in  how  CCTs  were  reported  and  (2)  the  older  practice  of 
collectively  analyzing  CCTs  calculated  from  spectral  radi¬ 
ances  and  irradiances. 

The  timing  and  duration  of  a  researcher’s  daylight 
measurements  clearly  affect  the  m(CCT)  calculated.  As 
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Henderson  points  out,  this  calculated  mode  should  coin¬ 
cide  with  that  actually  existing  at  the  observing  site.20 
In  other  words,  making  many  measurements  of  atypical 
spectra  will  bias  the  resulting  inverse  CCT  results.  Our 
long-term,  routine  acquisition  of  daylight  spectra  satisfies 
Henderson’s  recommendation,  and  thus  our  m(CCT)  re¬ 
sults  are  truly  representative  of  most  Granada  daylight. 

Although  Table  1  shows  that  clear  daylight  CCTs  have 
m(CCT)  ~  170-180MKT1,  the  CIE  in  fact  recommends 
using  the  D65  illuminant  with  CCT  =  6500  K  (or  154 
MKT1)  as  the  standard  illuminant  that  closely  approxi¬ 
mates  natural  daylight.  However,  the  timing  and  num¬ 
ber  of  our  clear-sky  daylight  measurements  bolster  our 
claim  that  our  m(CCT)  of  175-180  MKT1  is  indeed  the 
correct  one  for  our  site  rather  than  the  much  smaller 
m(CCT)  implied  by  the  CIE  standard. 

To  analyze  how  weather  affects  daylight  chromaticities 
in  Granada,  we  sorted  our  2600  spectra  into  three  bins 
that  represent  the  sky  states  listed  in  Table  2.  As  it  in¬ 
dicates,  clear  skies  occur  ~3.5  times  more  frequently  in 
Granada  than  do  overcasts.  Yet  perhaps  surprisingly, 
Fig.  4  shows  only  subtle  differences  in  inverse-CCT  fre¬ 
quency  distributions  for  these  two  extreme  sky  states. 


Nonetheless,  persistent  differences  do  exist:  m(CCT)  for 
overcasts  is  165-170  MK_1,  or  10  MK  1  smaller  (i.e., 
bluer)  than  that  for  clear  skies.  These  results  are  similar 
to  Dixon’s  in  the  southern  hemisphere.14  We  expect  this 
m(CCT)  difference  between  clear  and  overcast  skies:  The 
former  have  a  higher  ratio  of  yellow  sunlight  to  blue  sky¬ 
light  or  cloud  light  than  do  the  latter. 

Figure  5  shows  the  h0  dependence  of  CCT  for  clear  and 
overcast  skies.  The  mean  CCT  for  overcast  skies  is 
greater  than  that  for  clear  skies,  as  is  its  variability  (for 
h0  >  10°,  the  overcasts’  mean  CCT  exceeds  that  for  clear 
skies  at  the  5%  significance  level).  The  much  more  vari¬ 
able  nature  of  an  overcast’s  scattering  medium  means 
that  sudden,  large  changes  in  optical  depth  can  occur, 
with  resulting  fluctuations  in  chromaticity  and  irradiance 
beneath  the  overcast.  Figure  5  also  shows  that  neither 
clear  nor  overcast  skies  have  truly  meaningful  average 
CCTs  at  small  h0  (i.e.,  near  and  during  twilight)  because 
their  CCT  variance  is  so  large  then.  Note  too  that  mini¬ 
mum  CCT  occurs  for  20°  ^  Hq  ^  35  because  direct  sun¬ 
light  most  reddens  horizontal  irradiances  then.  At  still 
smaller  h0 ,  both  the  cosine  law  and  rapidly  increasing  at¬ 
mospheric  extinction  reduce  the  influence  on  horizontal 


Table  1.  Most-Frequent  Inverse  CCT  Intervals  m(CCT)  as  Measured  by  Various  Researchers 


Paper 

Sky  State 

m(CCT)  (in  MKT1) 

Henderson  and 

Hodgkiss1 

North  sky 

Total  sky,  no  sun 

Total  sky,  with  sun 

165-180 

145-150 

170-175 

Collins6 

Skylight  from  four  quadrants,  total  sky  with  sun 

165-170 

Winch  et  al. 1 

South  sky 

Total  sky,  no  sun 

Total  sky,  with  sun 

115-120,  160-165 
130-135 

170-175 

Ando  et  al.  (see 

Henderson20) 

North  sky,  clear  sky 

North  sky,  cloudy,  rain  or  snow 

95-100,  120-125 
160-165 

Tarrant11 

North  sky  and  south  sky 

150-160 

Dixon14 

Sunlight  plus  skylight: 

Partial  cloud  with  direct  sun 

Partial  cloud  with  obscured  sun 

Overcast  sky 

Clear  sky 

177.5- 182.5 

162.5- 167.5 

167.5- 172.5 

177.5- 182.5 

Present  work 

Hemispheric  daylight  (sunlight  when  present): 

Clear  sky 

Overcast  sky 

Few  clouds  to  mostly  cloudy  (i.e.,  all  others) 

175-180 

165-170 

175-180 

Table  2. 

Classification  of  our  Granada  Daylight  Spectra  by  Sky  State 

Sky  state 

Number 

of  Spectra  CCT  Range 

Inverse  CCT  Range 

Completely  clear 

Completely  overcast 

Few  clouds  to  mostly  cloudy 

757  5373  K-32,754  K 

216  3758  K-34,572  K 

1627  4864  K-26,728  K 

30.5  MK^-ise.i  Mr1 
28.9  MK_1-266.1  MKT1 
37.4  MK_1-205.6  MKT1 
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inverse  CCT  (MK1) 

Fig.  4.  Histogram  of  inverse  CCT  (in  MKT1)  for  Granada  daylight  measured  under  clear  and  overcast  skies. 


E(\  )  of  the  direct  beam’s  ever-reddening  spectra.  At  the  17000  y— - i - r 


largest  h0,  the  difference  between  maximum  and  mini-  isooo- 

mum  CCTs  is  only —200  K.  Still,  this  difference  exceeds  15000- 

1  JND.  14000- 

Finally,  we  looked  across  an  entire  year  for  possible  13000- 

seasonal  variations  in  Granada  daylight  chromaticities. 

^  12000  - 

While  a  1-year  survey  cannot  be  definitive,  we  found  no  ^ 
significant  shifts  either  in  our  locus  of  chromaticities  or  in  11000 

the  associated  histograms  of  inverse  CCT.  o  10000 ~ 


9000- 


■  clear  sky 
0  overcast  sky 


C.  Photometric  Characteristics 

In  many  applications,  we  need  photometric  rather  than 
radiometric  measures  of  daylight.  Because  atmospheric 
scattering  and  absorption  are  usually  spectrally  selective, 
ideally  one  calculates  illuminance  from  irradiance  rigor¬ 
ously  by  numerically  integrating  the  spectral  convolution 
of  E(\)  and  the  photopic  efficiency  function  Vp(\).32  A 
less  involved  (and  less  robust)  technique  is  to  use  a  single 
number  that  relates  illuminance  to  integrated  irradiance. 
This  number,  the  luminous  efficacy  factor,  is  the  ratio  of 
corresponding  photometric  and  radiometric  quantities  for 
a  particular  spectrum  (here,  the  ratio  of  illuminance  to  ir¬ 
radiance  integrated  over  the  solar  spectrum).  Such  a  fac¬ 
tor  lets  us  convert,  say,  a  pyranometer’s  irradiances  to  il¬ 
luminances. 

Given  our  spectroradiometer’s  spectral  range,  we  inte¬ 
grate  E(X)  from  300  to  1100  nm  when  calculating  lumi¬ 
nous  efficacy  factors.  Restricting  this  integration  to  the 
visible  (380-780  nm)  will  yield  significantly  larger  fac¬ 
tors,  since  by  definition  Vp(\)  =  0  outside  the  visible. 
Note  too  that  the  spectrum-dependent  luminous  efficacy 
factor  differs  from  the  constant  value  of  the  maximum 
photopic  luminous  efficacy,  Km  =  683  lumen/W. 32  As 
convenient  as  luminous  efficacy  factors  are,  using  a  vari¬ 
able  radiometric-to-photometric  conversion  clearly  poses 
potential  pitfalls.  Perez  et  al.  note  that  the  relationship 
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-5 


=5,  ±  =*=  ± 


“1 - 1 - 1 - 1 - 1 - 1 - r 

15  20  25  30  35  40  45 


“1 - r 

50  55 


n - r 

60  65 


hg  solar  elevation  (°) 


Fig.  5.  CCT  versus  solar  elevation  h0  at  5°  intervals  for  clear 
and  overcast  skies.  Squares  and  circles  indicate  mean  CCTs; 
each  error  bar  spans  2  standard  deviations. 


between  analogous  radiometric  and  photometric  quanti¬ 
ties  “still  leaves  many  questions  unanswered,  causing  po¬ 
tentially  large  errors  when  converting  one  type  of  quan¬ 
tity  into  the  other...”  (Ref.  33,  p.  33). 

Researchers  often  use  broadband  instruments  instead 
of  spectroradiometers  to  measure  outdoor 
illumination. 33-37  Broadband  devices  may  have  filters 
whose  spectral  transmissivity,  when  convolved  optically 
with  the  detector’s  spectral  sensitivity,  yields  a  net  spec¬ 
tral  response  that  closely  resembles  Vp(\).  However, 
because  no  broadband  filter-detector  combination 
matches  Vp(\)  exactly,  it  is  useful  to  develop  a  database 
of  daylight  illuminances  that  are  rigorously  calculated 
under  many  different  conditions.  For  a  given  daylight 
phase,  this  database  lets  a  researcher  accurately  convert 
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broadband  irradiance  to  illuminance  or  correct  a  photom¬ 
eter’s  fixed  illuminance  conversion. 

Figure  6  shows  such  a  database  of  luminous  efficacy 


Fig.  6.  Photopic  luminous  efficacies  [based  on  300- 1100  nm 
E  ( A )  integrals]  as  a  function  of  h0  for  our  2600  Granada  daylight 
spectra. 
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inverse  CCT  (MK1) 

Fig.  7.  Ratios  of  ultraviolet  (300-380  nm)  to  visible  (380-780 
nm)  irradiance  as  a  function  of  inverse  CCT  (in  MKT1).  The 
solid  curve  shows  these  ratios  for  Planckian  radiators  with  cor¬ 
responding  color  temperatures. 


Fig.  8.  Ratios  of  near-infrared  (780-1100  nm)  to  visible  (380- 
780  nm)  irradiance  as  a  function  of  inverse  CCT  (in  MKT1).  The 
solid  curve  shows  these  ratios  for  Planckian  radiators  with  cor¬ 
responding  color  temperatures. 


factors  calculated  from  our  Granada  daylight  spectra. 
These  factors  are  functions  of  h0  because  E(  \  )’s  shape 
changes  with  h0,  meaning  that  even  identical  integrated 
irradiances  may  produce  different  luminous  efficacies. 
For  h0>  20°,  luminous  efficacy  is  nearly  constant  at 
~140  lumen/W.  However,  it  increases  with  decreasing 
h0,  reaching  maximum  values  (and  maximum  variabil¬ 
ity)  when  the  sun  is  near  the  horizon  (h0  ~  5°)  and  the 
shapes  of  E(\)  and  Vp(\)  are  most  similar.  Luminous 
efficacy  plummets  after  sunset  as  twilight  E(X)  get  pro¬ 
gressively  bluer  [i.e.,  they  are  poorer  spectral  matches  to 
Vp(\)]. 

D.  UV  and  IR  Content 

Because  daylight’s  ultraviolet  component  can  profoundly 
affect  the  human  body,  being  able  to  infer  it  accurately 
from  visible-wavelength  irradiances  is  quite  useful. 
However,  Fig.  7  shows  that  such  inferences  are  problem¬ 
atic  at  best.  In  it,  for  each  Granada  daylight  spectrum 
we  calculate  the  ratio  of  ultraviolet  (300-380  nm)  to  vis¬ 
ible  (VIS)  (380-780  nm)  irradiances  and  then  plot  these 
UV/VIS  ratios  versus  inverse  CCT  for  the  visible  E(\). 
We  also  include  in  Fig.  7  as  a  reference  locus  the  UV/VIS 
ratios  for  Planckian  radiators  that  have  the  same  range  of 
color  temperatures  as  our  measured  E(  \ ) .  At  all  CCTs 
our  Granada  ratios  are  both  smaller  and  more  weakly  de¬ 
pendent  on  inverse  CCT  than  those  reported  by  Hender¬ 
son  and  Hodgkiss1  and  Winch  et  al.7  Yet  our  UV/VIS  ra¬ 
tios  do  resemble  those  of  Judd  et  al.2  at  low,  although  not 
high,  inverse  CCTs. 

Near-infrared  daylight  has  attracted  less  attention 
from  researchers  than  the  UV  (in  fact,  apparently  only 
Dixon  discusses  it).14  Historically  this  neglect  stems 
from  the  limited  (or  nonexistent)  IR  irradiances  in  most 
daylight  measurements.  To  redress  this  shortcoming,  we 
calculate  the  ratio  of  infrared  (780-1100  nm)  to  visible 
(380-780  nm)  irradiances  for  all  2600  of  our  Granada 
daylight  spectra.  As  shown  in  Figure  8,  these  IR/VIS  ra¬ 
tios’  correlation  with  inverse  CCT  is  much  larger  than 
that  for  the  UV/VIS  ratios.  In  fact,  as  Dixon  points  out14 
and  Fig.  8  reiterates,  the  IR/VIS  ratios’  relationship  to  in¬ 
verse  CCT  strongly  resembles  that  for  Planckian  radia¬ 
tors. 

Although  the  UV/VIS  and  IR/VIS  ratios  in  Figs.  7  and  8 
differ  from  some  reported  elsewhere,1’2’7’14  this  does  not 
necessarily  mean  that  the  UV  or  IR  content  of  Granada 
daylight  is  unusual.  Our  UV/VIS  and  IR/VIS  ratios  tell 
us  only  about  the  relative  balance  between  Granada’s  vis¬ 
ible  and  invisible  daylight  components,  and  this  balance 
depends  strongly  on  the  quantity  measured  (daylight  or 
skylight)  and  on  prevailing  atmospheric  conditions.  For 
example,  compared  with  skylight  spectra,  our  daylight 
spectra  usually  include  direct  sunlight  and  so  will  have 
proportionally  less  UV  energy.  Conversely,  our  spectra’s 
relative  IR  content  will  be  enhanced  whenever  they  in¬ 
clude  direct  sunlight. 

4.  SPECTRAL  ANALYSIS 

Color  analysis  and  synthesis  often  use  linear  models  to 
represent  illuminants.  When  paired  with  similar  reflec- 
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tance  models,  these  permit  the  design  of  efficient 
artificial-vision  algorithms.21,22  In  order  to  achieve  such 
efficiencies,  we  need  to  determine  the  smallest  dimension 
for  our  linear  daylight  models  that  lets  us  reconstruct  the 
original  E(\)  to  within  some  specified  accuracy.  To  de¬ 
termine  this  minimum  dimensionality,  we  perform  a 


Table  3.  Spectral  Regions  of  Our  2600  Granada 
Daylight  Spectra" 


Spectral  Region 

Spectral  Irradiances 

(nm) 

per  Region 

300-380 

17 

780-1100 

65 

380-780 

81 

300-1100 

161 

“Shown  are  the  number  of  spectral  irradiances  in  each  region  for  our 
spectroradiometer’s  5-nm  resolution. 


principal-value  decomposition  of  our  2600  Granada  day¬ 
light  spectra. 

We  begin  by  calculating  how  many  5-nm-wide  spectral 
bins  exist  in  each  of  the  spectral  regions  of  interest  for 
analyzing  daylight:  the  UV,  visible,  IR,  and  their  combi¬ 
nation  (Table  3).  Given  the  variability  of  daylight  in 
these  spectral  regions,  5-nm  spectral  resolution  yields 
more-reliable  reconstructions  than  does  10-nm 
resolution.24  Next  we  determine  the  eigenvectors  needed 
to  reconstruct  a  given  set  of  spectra  as  described  by 
Parkinnen  et  al.38  It  Ee(\)  is  a  measured  daylight  spec¬ 
trum,  then  we  can  approximate  it  by  using  p  eigenvectors 
in 

p 

Er(\)  =  2  (EeW |Vi(\)>Vi(k),  (2) 

i-l 

where  Er(K)  is  the  reconstructed  spectrum,  Vt( X)  is  the 
ith  eigenvector,  and  (|)  denotes  the  inner  product.  To  re- 


Fig.  9.  (a)  Spectral  distribution  from  300-1100  nm  of  eigenvectors  VpX),  V2(X),  and  V  3{\)  for  our  Granada  daylight  basis  set.  Solid 
curve,  i  =  1  eigenvector  VjfX)  [the  mean  2J(X)];  dashed  curve,  i  =  2  eigenvector  V2(X);  dotted  curve,  i  =  3  eigenvector  V3(X).  (b) 

Spectral  distribution  from  300  to  1100  nm  of  eigenvectors  V4(X),  W5( X. ) ,  and  P6(X)  for  our  Granada  daylight  basis  set.  Solid  curve,  i 
=  4  eigenvector  P4(X);  dashed  curve,  i  =  5  eigenvector  P5(X);  dotted  curve,  i  =  6  eigenvector  P6(X). 
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construct  a  curve  of  N  points  exactly,  in  principle  we  need 
N  eigenvectors.  In  practice,  however,  correlations  among 
our  2600  daylight  spectra  mean  that  we  can  set  p  <  N 
without  losing  any  meaningful  spectral  information. 

In  order  to  reconstruct  SPDs  with  any  desired  accu¬ 
racy,  we  must  adopt  the  approach  of  Eq.  (2)  rather  than 
use  the  approximate  CIE  technique.15  Only  the  former 
lets  us  adequately  analyze  the  dimensionality  p  needed  to 
achieve  a  given  accuracy  of  reconstruction.  To  test  our 
reconstructed  daylight  spectra,  we  once  again  use  a 
goodness-of-fit  coefficient17,23,24  (GFC)  that  is  based  on  the 
inequality  of  Schwartz.  Here  this  coefficient  is 

=  I2/ 

|S,[jE£(Aj)]2|1/2|2,[^(^)]2|1/2' 

GFC  is  the  multiple  correlation  coefficient  R,  the  square 
root  of  E^Ol’s  spectral  variance  with  respect  to  the  origi¬ 
nal  Ee(\).  GFC  ranges  from  0  to  1,  where  1  indicates  a 
perfect  reconstruction.  We  found  that  colorimetrically 
accurate  ER(\)  require  GFC  3=  0.995;  what  we  call  “good” 
spectral  fit  requires  GFC  3=  0.999,  and  GFC  3=  0.9999  is 
necessary  for  “excellent”  or  almost-exact  fit.  Examples  of 
Ee(\)  with  different  GFCs  are  found  in  Refs.  17  and  24. 

In  Fig.  9  we  plot  the  first  six  Granada  daylight  V;(X) 
for  300-1100  nm.  The  i  —  1  eigenvector  Vj  ( X)  corre¬ 
sponds  to  our  database’s  mean  £(X),  and  its  spectral  de¬ 
tails  include  only  a  few  small  absorption  bands.  As  Judd 
et  al.2  and  Sastri  and  Das9  note,  the  V2(X)  eigenvector  ac¬ 
counts  for  most  yellow-blue  chromaticity  variations. 
Like  the  mean  E{\),  V2(X)  also  varies  smoothly,  increas¬ 
ing  from  negative  to  positive  values  near  X.  =  550  nm. 
Starkly  different  is  the  V3(X)  eigenvector,  many  of  whose 
sharp  peaks  derive  from  atmospheric  absorption  bands. 
Here  we  part  company  with  Judd  et  al.  and  Sastri  and 


Das:  Unlike  them,  we  find  no  purple-green  shift  associ¬ 
ated  with  V3(X).  As  Fig.  10  shows,  V3( X )  actually  yields 
chromaticity  shifts  nearly  parallel  to  those  of  V2(X).  To 
determine  the  chromaticity  shifts  plotted  in  Fig.  10,  we 
calculate  the  chromaticity  pairs  that  result  when  we  add 
and  subtract  either  V2(X)  or  V3(X)  to  and  from  Fj(X). 

Although  our  collection  of  daylight  spectra  is  the  larg¬ 
est  measured  to  date,  we  require  only  a  few  V,(X)  to  ac¬ 
count  for  nearly  all  its  variance  in  the  visible.  Unlike 
earlier  researchers,4,9  we  find  that  the  first  five  F;(X)  ac¬ 
count  for  99.991%  of  the  observed  variance  between  380 
and  780  nm.  This  good  performance  in  the  visible  is  not 


Fig.  10.  Colorimetric  characteristics  of  the  eigenvectors  shown 
in  Fig.  10(a),  with  the  chromaticity  of  the  mean  vector  TJX) 
marked  “+.”  For  V2(X )  and  V3(X)  we  plot  the  chromaticity  shift 
that  each  accounts  for  (lines  with  arrows).  We  also  include  the 
deuteranopic  and  tritanopic  confusion  lines  (skewed  grid  lines) 
as  well  as  the  Planckian  (dashed  curve). 


Table  4.  Percentage  of  Reconstructed  Spectra  ER(\.)  That  Exceed  Our  Three  Target  GFC  Values  As  a 
Function  of  Spectral  Region  and  Number  of  Eigenvectors  p  Used  in  Eq.  (2)  (Here,  2  =£  p  =£  8)“ 


Spectral 

Region  (nm) 
[Number  of 

E(\) 

per  spectrum] 

p  Eigenvectors 

2 

3 

4 

5 

6 

7 

8 

300-380 

99.7 

99.9 

100 

(17) 

95.9 

98.0 

98.5 

98.7 

98.8 

99.3 

99.4 

37.4 

78.6 

93.4 

93.6 

94.5 

94.9 

95.3 

780-1100 

90.0 

99.8 

100 

(65) 

54.6 

98.0 

99.3 

99.5 

99.8 

99.8 

99.8 

8.4 

43.0 

80.5 

93.0 

95.3 

97.2 

97.4 

380-780 

95.2 

99.8 

100 

(81) 

78.0 

97.4 

99.5 

99.9 

100 

3.8 

29.3 

61.5 

91.5 

97.2 

99.1 

99.4 

300-1100 

92.5 

98.3 

99.5 

99.9 

100 

(161) 

37.9 

70.7 

88.5 

97.1 

99.6 

99.8 

99.8 

0.3 

4.5 

17.5 

48.7 

70.7 

88.0 

92.3 

Within  each  cell,  rows  from  top  to  bottom  are  for  target  GFC  =  0.995,  0.999,  and  0.9999. 
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surprising  given  the  Granada  chromaticities’  close  prox¬ 
imity  to  their  daylight  locus  (Fig.  2).  Elsewhere  in  the 
spectrum,  however,  the  E(X)  correlation  is  not  as  high. 
For  example,  the  correlation  is  lower  from  300  to  380  nm 
(which  includes  UVA  and  UVB),  although  reduced  accu¬ 
racy  of  our  spectroradiometer  here  may  increase  variance. 
As  in  the  visible,  from  300  to  380  nm  the  first  five  V^X) 
account  for  99.990%  of  the  total  observed  variance,  even 
though  the  visible  spans  a  wavelength  region  five  times 
larger  (and  thus  p  —  5  might  be  expected  to  account  for 
less  variance  there).  In  the  LI-1800’s  320-nm  wide  IR 
band,  the  first  five  V^X)  account  for  about  as  much  of  the 
total  variance  as  they  do  in  the  UV. 

Next  we  use  the  CIE  method15  to  estimate  daylight 
Ee(X)  from  their  associated  chromaticities.  However,  we 
update  this  method  by  using  the  Granada  V,-(\)  shown  in 
Fig.  9(a)  and  Eq.  (l)’s  Granada  daylight  locus.  Because 
we  reconstruct  ER(X)  using  p  =  3,  the  mean  GFC  ob¬ 
tained  for  our  2600  ER(X)  is  only  0.99543.  Moreover,  we 
note  that  merely  33%  of  these  ER(X)  have  GFC  >  0.999 
and  only  4.6%  have  GFC  >  0.9999. 

Using  2  =s  p  =£  8  in  Eq.  (2),  we  calculated  all  2600 
Er( X)  for  Table  3’s  four  spectral  ranges.  Table  4  and 
Fig.  11  indicate  the  accuracies  possible  for  these  upper 
limits.  Not  surprisingly,  mean  GFC  increases  with  p,  al¬ 
though  the  increase  is  not  uniform  in  all  spectral  regions. 
In  particular,  with  p  =  2  the  UV  yields  the  highest  mean 
GFC  compared  with  the  remaining  spectral  regions.  By 
contrast,  using  p  >  5  scarcely  increases  mean  GFC  in  the 
UV  from  its  p  =  4  values. 

Figure  11  and  Table  4  indicate  that  the  visible  pro¬ 
duces  the  best  mean  GFC  and  that  the  greatest  variabil¬ 
ity  (smallest  correlation)  exists  in  the  UV  and  IR.  For 
the  mean  GFC  in  all  four  spectral  ranges  to  be  considered 
“good,”  we  require  p  s*  3,  and  for  it  to  be  “almost  exact” 
we  need  p  s®  6.  However,  we  should  not  choose  p  solely 
on  the  basis  of  mean  GFC,  independent  of  spectral  region. 
Instead,  we  need  to  see  how  few  Vt(X)  are  required  for  us 
to  meet  our  target  GFCs  in  each  spectral  region  (Table  4). 
For  95%  of  Er(X)  to  have  at  least  “good”  quality,  the  mini¬ 
mum  p  varies  from  2  in  the  narrowest  region  (the  UV)  to 
5  in  the  widest  (300-1100  nm).  For  95%  of  ER(X )  to  be 
“almost  exact,”  p  may  need  to  be  as  high  as  10  from  300  to 
1100  nm,  8  in  the  UV,  or  6  in  the  visible  and  IR. 

The  CIE  recommends  using  only  three  basis  vectors  (al¬ 
though  with  2  degrees  of  freedom  each)  to  achieve  a 
univocal  relationship  between  chromaticities  and  recon¬ 
structed  E{  X ) .  Yet  if  we  want  a  spectrally  “almost-exact” 
Er(X)  at  visible  wavelengths,  we  must  use  p  s®  6  in  Eq. 
(2).  Adding  either  the  UV  or  the  IR  region  to  the  visible 
requires  p  >  7  (as  Slater  and  Healey  note).16  If  we  want 
highly  accurate  E R(X)  from  300  to  1100  nm,  we  must  use 

p  =  10. 

Although  increasing  p  always  improves  ER(X)’s  aver¬ 
age  spectral  quality,  note  that  this  extra  computational 
burden  soon  produces  only  marginal  colorimetric  im¬ 
provements  in  Er,(X).  Because  the  relationship  between 
Er OO’s  spectral  and  colorimetric  accuracy  is  not  simple, 
we  perform  yet  another  test.  In  order  to  evaluate  the 
size  of  color  differences,  we  convert  our  colorimetric  data 
to  the  CIELUV  uniform  color  space.39  In  calculating  the 
CIELUV  color  difference  AL*(J  we  use40 


number  of  eigenvectors  p 

Fig.  11.  Mean  GFC  for  our  2600  Granada  daylight  spectra,  with 
use  of  2  p  10  eigenvectors  in  Eq.  (2)  for  the  indicated  spec¬ 
tral  region.  We  label  in  parentheses  the  number  of  spectral  ir- 
radiances  within  each  region  for  5-nm  resolution. 


Fig.  12.  Mean  CIELUV  color  difference  tXE*v  (Eq.  4)  as  a  func¬ 
tion  of  mean  GFC.  Numbers  next  to  the  squares  indicate  the 
number  of  eigenvectors  p  needed  in  Eq.  (2)  to  obtain  the  given 
GFC.  Each  error  bar  spans  2  standard  deviations. 

A  E*v  =  [(AL*)2  +  (Air*)2  +  (Au*)2]1/2,  (4) 

where  L*,  u*,  and  u*  are  the  CIELUV  space’s  orthogonal 
coordinates,  and  AL*,  Am*,  A u*  are  the  corresponding 
differences  between  coordinates  of  the  original  and  the  re¬ 
constructed  colors. 

Figure  12  shows  how  the  mean  CIELUV  color  differ¬ 
ence  between  Ee(X)  and  ER(X)  changes  as  p  increases 
from  2  to  10  (the  mean  here  is  over  all  2600  daylight  spec¬ 
tra).  As  p  increases,  the  mean  colorimetric  distance 
A E*v  decreases  sharply,  especially  for  p  <  5.  Three  or 
five  A E*v  units  are  often  taken  to  be  1  JND  in  technical 
and  industrial  applications.  Thus  using  p  =  3  in  Eq.  (2) 
means  that  both  the  mean  A E*v  and  its  standard  devia¬ 
tion  will  satisfy  this  practical  JND  criterion.  So  from  a 
colorimetric  (rather  than  a  spectral)  standpoint,  the  CIE’s 
recommendation  of  three  basis  vectors  is  adequate.  Note 
that  we  reached  this  conclusion  by  having  the  original 
Ee (X)  available  to  use  in  Eq.  (2).  Yet  when  we  actually 
implement  color  identification  and  recognition  algo¬ 
rithms,  we  will  lack  the  original  Ee{X)  and  instead  must 
rely  on  sensors  with  poorer  spectral  resolution  and  per¬ 
haps  additional  noise.  Thus  for  truly  reliable  results, 
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prudence  suggests  that  we  use  more  than  the  absolute 
minimum  number  of  V,-(\)  in  Eq.  (2),  with  p  —  4  striking 
a  good  compromise  between  efficiency  and  error. 

5.  CONCLUSIONS 

Our  suite  of  horizontal  daylight  E(\)  provides  an  excep¬ 
tionally  broad  range  of  natural-light  chromaticities  and 
CCTs.  The  differences  between  the  CIE  and  the 
Granada  daylight  loci  are  most  evident  at  high  CCTs 
(>9000  K),  where  the  Granada  locus  is  greener.  If  we 
had  to  recommend  one  CCT  as  best  characterizing 
Granada  daylight,  our  analysis  suggests  5700  K.  Our 
measurements  also  indicate  that  Granada’s  most- 
frequent  inverse  CCT  [m(CCT)]  under  clear  skies  is  175- 
180  MKT1  (or  5555-5715  K),  a  value  that  decreases  to 
165-170  MK -1  (or  5882-6060  K)  under  overcasts. 
Granada’s  UV/VIS  ratios  are  less  dependent  on  CCT  than 
any  reported  previously;  the  analogous  IR/VIS  ratios 
nearly  duplicate  those  of  Planckian  radiators  with  the 
same  color  temperatures. 

A  principal-components  analysis  of  our  2600  Granada 
daylight  spectra  indicates  that  the  i  —  2  eigenvector  ac¬ 
counts  for  their  yellow-blue  chromaticity  variations  and 
that  the  i  =  3  eigenvector  does  not  yield  any  purple- 
green  variation,  although  it  does  depend  on  atmospheric 
absorption  bands.  Spectrally  speaking,  our  work  shows 
that  we  need p  >  5  for  accurate  E  ft( k )  in  the  visible,  even 
though  CIE  recommendations  for  reconstructing  daylight 
SPDs  from  chromaticities  in  effect  call  only  for  p  —  3.  If 
we  include  either  the  UV  or  the  IR  in  our  ER(\),  the  mini¬ 
mum  p  increases  to  7;  if  we  want  to  construct  the  entire 
300-1100  nm  spectrum,  we  need  10  vectors.  As  the  CIE 
recommends,  p  —  3  will  produce  ErCK)  that  are  colori- 
metrically  indistinguishable  from  the  corresponding 
Ee(X)  for  most  observers,  although  in  machine-vision  ap¬ 
plications  a  modest  increase  to  p  —  4  seems  advisable. 

We  would  not  be  surprised  if  subsequent  daylight  mea¬ 
surement  campaigns  in  other  countries  produce  colori¬ 
metric  and  spectral  analyses  that  differ  perceptibly  from 
ours.  However,  our  earlier  work  suggests  that  daylight 
spectra  do  share  some  fundamental,  quantifiable  similari¬ 
ties  worldwide.24  Thus  we  view  this  paper  as  a  vital  first 
step  in  determining  whether  and  to  what  extent  site- 
dependent  differences  are  significant  in  daylight.  More 
immediately,  our  analysis  of  the  largest  extant  set  of 
hemispheric  daylight  spectra  offers  useful  baseline  infor¬ 
mation  for  any  site  that  lacks  comparable  data.  As  is 
true  elsewhere  in  atmospheric  optics,  revisiting  the  well- 
trod  ground  of  daylight  spectra  still  affords  fresh  spectral 
and  colorimetric  insights.  In  addition  to  providing  a 
large,  valuable  library  of  such  data,  we  believe  that  our 
work  suggests  new  ways  of  understanding  and  solving 
some  quite  practical  problems  associated  with  outdoor  vi¬ 
sion,  whether  machine  or  human. 
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